J/ip polarization in photoproduction up-to the next-to-leading order of QCD 
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We investigate the J/ if) polarization in photoproduction at the Hadron Electron Ring Accelera- 
tor(HERA) up to the next-to-leading order of QCD. The results show that the transverse momentum 
Pt and energy fraction z distributions of J/'tp production do not agree with the observed ones very 
well. The theoretical uncertainties for the z distributions of the .J/^j) polarization parameters with 
respect to various choices of the renormalization and factorization scales are too large to give an 
accurate prediction. The uncertainties for the pt distributions of these parameters are small when 
Pt > 3GeV and the obtained pt distributions can not describe the experimental data even in this 
region. 
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The heavy-quarkonium systems (cc) and (66), espe- 
cially J ftp and T, being flavor hidden, have received 
much attention from both experimental and theoretical 
sides, since they were discovered. Their heavy masses set 
a "large" scale that makes the effective field theory, non- 
relativistic quantum chromodynamics (NRQCD)[l| ap- 
plicable for their production and decay processes. In 
terms of factorization, NRQCD manages the expansion 
on the strong coupling constant as in the hard part and 
the V (velocity between the heavy quark and antiquark in 
quarkonium) in the soft part properly. Indeed, NRQCD 
achieves much success by taking high Fock state contribu- 
tions into account. It reconciles several discrepancies be- 
tween the theoretical predictions and experimental data, 
but there are still some problems. Of the problems, the 
discrepancies between the theoretical predictions at lead- 
ing order (LO) and experimental data on the polarization 
of J/ip and T hadroproduction are outstanding. Recent 
reviews on the situation can be found in Ref. [2|. 

Many works 0, i, i] indicate that the next-to-leading- 
order (NLO) QCD corrections under the NRQCD frame- 
work drastically change the features of LO theoretical 
predictions in the heavy-quarkonium production in vari- 
ous cases. Especially, it is reported in Ref.Q that with 
the color-singlet model (CSM) alone the NLO predic- 
tion on J/ip photoproduction at HERA can explain the 
experimental data well. Also with the CSM alone, the 
NLO results on the charmonium production moderate, 
even resolve, the discrepancies between theoretical pre- 
dictions and experimental measurements at B factories 
Recently, the NLO QCD corrections on the hadronic 
production of J/ip and T in the CSM have also been stud- 
ied by several groups @, 0, [1] and the results show that 
the Pt distribution of J/^ production is largely enhanced. 
Furthermore, the NLO QCD corrections on the hadronic 
production of J/V' in the color-octet mechanism (COM) 
have been completed in Ref. [9| and the results show that 
the Pt distribution of J/ip production is changed slightly, 
whereas the polarization of the produced quarkonium 
still is an open problem. Namely, the results at NLO 
in the CSM give a longitudinal polarization contrary to 



the transverse one as LO doespT^, and the NLO results in 
the COM give almost unchanged polarization as the LO 
ones [i]. To add the NLO CSM and NLO COM results 
together, the data for polarization of J/ip hadroproduc- 
tion at Tevatron cannot be described properly, although 
the data forp* distribution of hadroproduction can be fit 
quite well li- 
lt may provide ep and 7p processes to investigate the 
production of J/ip at HERA. The transverse momentum 
Pt distributions of J/ip production and polarization at LO 
were studied in the CSM a long time ago . 
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As mentioned above, it was reported in Ref.[3| that with 
the CSM alone the NLO results can give a proper descrip- 
tion for the photoproduction of J/tp at HERA, including 
the energy fraction distribution da/dz in the intermedi- 
ate z region and transverse momentum pt distribution, 
where the energy fraction z is defined by z = ^"^/^"^ ^''^ 

and pj/^, p-y and pp are the momenta of J/ 4', 7 and pro- 
ton, respectively. Then, to investigate the J/ip inelastic 
photoproduction, many works based on NRQCD at LO 
and NLO in the CSM and COM followed[i3]. In Ref. [3 
the photoproduction of polarized JM was studied at LO 
in the CSM and COM. In Ref. [3 the theoretical pre- 
dictions based on the kx factorization formula were also 
given at LO. However, the latest experimental results on 
the photoproduction of polarized J/V'tlZi] do not favor 
the LO predictions in the CSM and COM, even those ob- 
tained with the kt factorization formula. So, up to now 
there is no satisfactory explanation on the measurement 
results of photoproduction of polarized J/tp at HERA. 
In view of the proper description of the HERA data on 
Pt and z distributions by the CSM NLO prediction^, 
we focus the theoretical predictions on J/tp polarization 
with NLO QCD corrections in the CSM in this paper. 

To calculate the NLO QCD correction for the photo- 
production of polarized J/tp in the CSM, there are direct 
and resolved processes to be considered. The existing 
LO calculation [ll] shows that the contribution from re- 
solved processes is about 2 orders of magnitude smaller 
than that of the direct one for the pt distribution of J/tp, 
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and only in the lower z (z<0.2) region is the contribution 
from resolved processes comparable with that from the 
direct one. Therefore, with the cut condition z > 0.2 
for the experimental measurements [itI [18|. only the di- 
rect processes are investigated in this paper. To perform 
the lengthy analytic evaluation for all the processes, the 
computer program package, Feynman Diagram Calcula- 
tion (FDC), is used. FDC was developed and well tested 
for many years [l^, and recently the functions for ma- 
nipulating one-loop quantum corrections in analytic re- 
duction and numerical calculations were completed and 
tested in many aspects j20|. 

To calculate the J/ip photoproduction at NLO, the 
following relevant processes 



7 + 5^ J/ip + 9, 
1 + 9^ J/-0 + g + g, 
1 + 9 ^ J/ij + q + q, 
i + q{q) J/iP + 9 + q{q)- 



(1) 
(2) 

(3) 
(4) 



need to be considered. For process ([T]), there are 6 Feyn- 
man diagrams at LO and 111 at NLO, and there are 
the ultraviolet (UV), infrared(IR) and Coulomb singular- 
ities. By using dimensional regularization and introduc- 
ing a small relative velocity between the quark and an- 
tiquark, we can separate the singularities out. The UV 
and Coulomb singularities are absorbed into the redef- 
inition of mass, coupling constant, fields of quark and 
gluons and wave function of J/tp by renormalization, for 
which the same renormalization scheme as in R ef. lH is 
applied. By using the phase space shcing method [2l|. we 
separate out the IR singularities in the real processes ^ 
([3]) ([4]). Finally to add all the contributions from the real 
and virtual parts together, a finite result, which is free 
from all the singularities, is obtained. 

As the main purpose of this paper, we extract the in- 
formation of the J/f/' polarization through the angular 
distribution of the lepton in J/tp decay. As in 

Ref. . the decay angular distribution of the outgoing 
Z"*" can be parametrized in the J /ip rest frame as 



da 
dfldy 



(X 1 + \{y) cos^ 6 + fi{y) sin 29 cos (/) 
+ '^sm^ 9cos2(f>, 



(5) 



where y stands for a suitable variable (such as transverse 
momentum pt, energy fraction z, etc.). 9 and are the 
polar and azimuthal angles of the outgoing respec- 
tively. The polarization parameters. A, /i, and v, are 
related to the density matrix of J/ijj production as 



Ky)-- 
Ky) = 



daii/dy ~ dapo/dy 
daii/dy + daoo/dy' 

\/2Redaio / dy 
daii/dy + daoo/dy' 

2dai-i/dy 
daii/dy + daoo/dy' 



(6) 



Here daw /dy are the differential density matrix ele- 
ments and defined as 



daxx> 
dy 
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5{y~y{p,))M{X)M*{\'), 



(7) 



where M{\) is the matrix element of polarized J /ip pro- 
duction, A and A' stand for the polarization, F includes 
the flux factor and spin average factor, and pi is the mo- 
mentum of corresponding particles. In the calculation, 
the polarization of J /ip must be explicitly retained and 
the treatment for it is the same as in Ref. Q • Obviously 
these polarization parameters depend on the coordinate 
system choice. In the following, we calculate them in 
the target frame with the Z axis defined as the inverse 
direction of the initial proton and the polarization vec- 
tors of J/%p defined in the appendix of Ref. [l^. We 
also compute the polarization parameters of J /^ in the 
helicity-base frame where the Z axis is defined as the J /"ip 
fiight direction in the laboratory frame, and the angular 
distribution of is parametrized as 



da 



d cos 9dy 



oc 1 -I- a{y) cos^ 



(8) 



Here the parameter a is related to the polarized cross 
sections of J/?/' production by 



a{y) 



dar/dy - 2daL/dy 
dar/dy + 2daL/dy 



(9) 



where (Tt and ctl are the cross sections of transverse and 
longitudinal polarized J/tp respectively, a = —1 corre- 
sponds to fully longitudinal polarization and a = 1 to 
fully transverse polarization. 

To replace the polarization vectors for photon or glu- 
ons by their corresponding momentum in the numerical 
calculation, the gauge invariance is obviously observed. 
Since the two phase space cutoffs are chosen to handle 
the IR singularities of the real processes, we numerically 
check the independence of the results on the cutoffs. 

In the numerical calculation, we use q;=1/137, 
mc=1.5GeV, and Mj/^ = 2mc- The CteqGLl and 
Cteq6M[23 are used in the calculations at LO and NLO 
respectively, with the corresponding as running formula 
in the Cteq6 being used. In Ref. a fixed value of 
as is taken for numerical calculations. Whereas the 
running coupling constant is chosen in most literatures, 
here we use the running as in all our calculations. The 
renormalization scale fir and the factorization scale /z/ 

are set as fir = fJ-f = \/ (2mc)^ +Pt • The value for 
J/ip wave function at the origin is extracted from the 
leptonic decay width with the NLO formula, the value 
Fj/^^ee ==5.55keV, and as(Mj/^) =0.26. The typi- 
cal HERA photon-proton center-of-mass energy ,Js^ = 
lOOGeV is chosen. In addition, the experimental cut con- 
ditions 0.4 < z < 0.9 and pt>\ for the pt and z distri- 
butions on J /ip production are applied respectively, and 
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FIG. 1: The energy fraction z distribution of J/?/) photoproduction with Pt >lGeV (left panel). The transverse momentum pt 
distribution of J/?/) with 0.4 < z < 0.9 (right panel). The dotted lines in each one are the LO results. The upper bound of the 
shaded band is obtained with ~ — \ \/ (2mc)^ and mc=1.4GeV, and the lower one with fir ~ fJ-f = 2^ (2mc)^ + Pt 
and mc = 1.6GeV. The experimental data (filled circles with bars) are taken from Ref. [l^ . 
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FIG, 2: The energy fraction z distributions of polarization parameters A, v, and fi with pt > IGeV, and transverse momentum 
Pt distributions with 0,4 < z < 1. Dashed lines are the LO results and the dot-dashed lines are the NLO results. Solid lines 
present the results with fir = fJ-f = Srric. The experimental data (filled circles with bars) are taken from Ref. [l7| ]. 



0.4 < z < 1 and pt > 1 are applied for the J/ip polariza- 
tion. Furthermore, all the differential cross sections are 
calculated directly in numerical calculation by using the 
corresponding analytic phase space treatment just like 
that used in Ref, 

The final results are presented in the figures. From 
FiglU we can see that the results for pt and z distribu- 
tions of J/tp production at NLO moderate the discrep- 
ancies between theoretical predictions and experimental 
measurements substantially. The pt and z distributions 
of the polarization parameters in Eq,([5]) is presented in 
Fig, [2] In the presentation, we make a cutoff for the re- 
gion z < 0,2, because in this region the denominator of 
Eq, [S] crosses zero at a certain point so that the values of 
the parameters change dramatically in the neighborhood 
of the point. It means that the perturbation expansion 



is very bad and the obtained result cannot be trusted in 
this region. Comparing with the results at LO, the NLO 
QCD corrections greatly change the distributions. The 
parameters A obtained at LO and NLO are very different 
at the large and small z regions (there is a cross at about 
z « 0,67), The z distribution of parameter i' changes 
drastically in the small and intermediate regions, and fi 
also changes in the intermediate region of z. As for the pt 
distribution of the parameters, from the figure one can 
see that the NLO QCD correction changes that of the 
A parameter from positive values to negative ones and 
makes it tend to -0,8 as pt increases. For the v param- 
eter, there is little difference between the results of LO 
and NLO in the large pt region. The influence of the 
NLO correction on the pt distribution of fi is also quite 
large. For comparison, the available experimental data at 
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FIG. 3: The polarization parameter a distributions as func- 
tions of Pt and z in the helicity basis. Dashed hnes are the 
results at LO and the dot-dashed lines are the ones at NLO. 
The results with fi,, = fJ^f ~ 8mc are presented by the solid 
lines. 



HERA [iTi are plotted in the figures. It is clearly shown 
that both LO and NLO results do not fit the polariza- 
tion measurement, and the NLO results are even worse. 
In these figures, the z and pt distributions of A,^, and v 
with = fif = Snic at NLO are also presented. 

The Pt and z distributions of a in the helicity basis are 
presented in Fig. [3] We also make a cutoff for the region 
z < 0.2, due to the same reason as that in the case of 
the target frame. It is similar to the situation in the J/tp 
hadroproduction: NLO QCD corrections change the pt 
distribution of a from positive values to negative ones. 
The z distribution of a changes similar to that on the pt 
distribution except the endpoint region near z — 1. We 
also present the results with = = 8mc. 

In summary, we have investigated the photoproduction 
of J/ip at QCD NLO at HERA. The results show that 
the NLO QCD corrections cannot give a very good de- 

^ Private communication with P. Artoisenet and F. Maltoni 



scription on the pt and z distributions. Furthermore, for 
the J/tp polarization, the NLO QCD corrections make 
a drastic change on certain polarization parameters, but 
they cannot give a satisfied prescription for the available 
experimental data. In contrast to the LO calculation 
for polarization, the NLO results are even worse. There 
are two conclusions on the results at QCD NLO. One 
is that the theoretical uncertainties for the z distribu- 
tions of the J/ip polarization parameters A, /z, u, and a 
on different choices of the renormalization and factoriza- 
tion scales are too large to give reasonable predictions to 
compare with the experimental measurement. Another 
is that the theoretical uncertainties for the pt distribu- 
tions of A, fi, V, and a on different choices of the scales 
are small when pt > iGeV and the obtained pt distri- 
butions cannot describe the experimental data even just 
in this region. Therefore, there is still no satisfactory 
theoretical description even at the NLO level on the pt 
or z distributions of J/ if) polarization as well as the pt 
or z distribution of J/i/' productions at HERA. It will 
be interesting to know the situation by considering NLO 
color-octet contributions, or higher order QCD correc- 
tions in a future study. 

While this paper was being prepared, we were informed 
of the same process also being considered by Artoisenet 
et al. [23]. Comparing our results in Figs. 1 and 2 
with theirs, there are quantitative discrepancies between 
them. But our results on the z distributions of A and v 
with /i,. = = StUc at NLO are consistent with theirs 
with y/Sjp — lOOGeV and /i^ = /^/ = Sruc^. It can be 
inferred that the discrepancies mainly come from the dif- 
ferent choices of renormalization and factorization scales. 
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